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ABSTRACT 

Many of the results in modern astrophysics rest on the notion that the Initial Mass Function (IMF) is uni- 
versal. Our observations of a sample of H I selected galaxies in the light of Ha and the far-ultraviolet (FUV) 
challenge this result. The extinction corrected flux ratio Fu a /fFuv from these two tracers of star formation 
shows strong correlations with the surface-brightness in Ha and the R band: Low Surface Brightness (LSB) 
galaxies have lower Fhq.//fuv ratios compared to High Surface Brightness (HSB) galaxies as well as compared 
to expectations from equilibrium models of constant star formation rate (SFR) using commonly favored IMF 
parameters. Weaker but significant correlations of Fh q //fuv with luminosity, rotational velocity and dynami- 
cal mass are found as well as a systematic trend with morphology. The correlated variations of fka//FUV with 
other global parameters are thus part of the larger family of galaxy scaling relations. The fka/ZFUv correla- 
tions can not be due to residual extinction correction errors, while systematic variations in the star formation 
history can not explain the trends with both Ha and R surface brightness nor with other global properties. The 
possibility that LSB galaxies have a higher escape fraction of ionizing photons seems inconsistent with their 
high gas fraction, and observations of color-magnitude diagrams of a few systems which indicate a real deficit 
of O stars. The most plausible explanation for the correlations is the systematic variations of the upper mass 
limit M. u and/or the slope 7 which define the upper end of the IMF. We outline a scenario of pressure driving 
the correlations by setting the efficiency of the formation of the dense star clusters where the highest mass stars 
preferentially form. Our results imply that the star formation rate measured in a galaxy is highly sensitive to 
the tracer used in the measurement. A non-universal IMF would also call into question interpretation of metal 
abundance patterns in dwarf galaxies as well star formation histories derived from color magnitude diagrams. 
Subject headings: galaxies: stellar content - stars: luminosity function, mass function - stars: formation - 
galaxies: ISM 
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1. INTRODUCTION 

A key concept in our understanding of the evolution of 
galaxies is the Initial Mass Function (IMF), which gives the 
statistical distribution of masses of stars that form in a sin- 
gle event. Typically the IMF is parametrized as a power law, 
or a series of broken power-laws, in stellar mass, where 
the key parameters are the power-law index 7 and the lower 
and upper mass limits to the stars formed Mi and J\A U respec- 
tively. The IMF is crucial in interpreting the nature of galaxies 
because the light we observe is dominated by the highest mass 
stars while the total mass in stars is dominated by the lower 
mass stars. If we know the IMF of a stellar population, and 
its age, then we can estimate its mass from photometry. Simi- 
larly when looking at tracers of star formation, such as Ho or 
ultraviolet (UV) emission, the IMF allows us to convert from 
luminosity to a star formation rate (SFR). 

The crucial assumption un derpinning mu ch of the use of 
the IMF is that it is universal (Gilmore 2001); it does not vary 
within galaxies or between galaxies. This certai nly seems to 
be the case for star clusters (Kroupa 200Tll2002h. Since most 
star formation occurs in star clusters dLada & La da 2003) a 
universal IMF seems plausible, and is well accepted in the 
literature. 

In recent years there has been accumu lating evidence tha t 
the IMF is not universal (see review of Elmegreen 2008a). 
This sets the stage for our study in which we challenge the 
notion that the upper end of the IMF is uniform using mea- 
surements of the integrated Ha and UV emission in a sam- 
ple of H I selected galaxies from the Survey of Ionization in 
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Neutral Gas Galaxies (SINGG), and the Survey of Ultravi- 
olet emission in Neutral Gas Galaxies (SUNGG). Ha emis- 
sion traces the presence of ionizing O stars, which have initial 
masses > 20.M©, while UV emission traces both O and 
B stars with Ai+ > 3A4q, giving the leverage in mass needed 
to probe the IMF. 

The H I selection performed by SINGG and SUNGG re- 
sults in a sample of star forming galaxies spanning the full 
range of morphologies se en in star forming g alaxies with very 
few Ha non-detections (Meu rer et al.l 12006. hereafter M06) 
and can fully account for the cosmic den sity of star formation 
in the local universe (lHanish et al.l l200ft hereafter H06). By 
using this sample we avoid the issue of selection effects bias- 
ing towards particular extreme star formation histories (SFH) 
such as starburst galaxies. While the SFH of the universe as 
whole appears to have undergone a dramatic factor of 10 de- 
cline since z ~ 1 (~ 8 Gyr ago for Hq = 70kms~ 1 Mpc~ 1 which 
we adopt here) this is still long compared to the lifetimes of 
the O (< 7 Myr) and B (< 320 Myr) stars probed in this study, 
so the SFH should be effectively constant for the sample as a 
whole. 

There are a few other advantages of our approach. We 
use int egrated fluxes, so in the parlance of Weidne r" & Kroupal 
J2005I) we are dealing with the Integrated Galaxial IMF 
(IGIMF). Unlike the case of color-magnitude diagram anal- 
ysis, crowding effects are irrelevant as are concerns as to the 
birth place of stars and their eventual dispersal in the field. 
Stochastic effects are minimized by integrating over whole 
galaxies, and the results are more comparable to studies of 
distant galaxies for which little more than total fluxes are 
available. 

This paper is organized as follows. In Sec. [2] we discuss the 
sample selection, the data used, the measurements of these 
data which we make, and the necessary corrections to the 
measurements. In Sec.[3]we demonstrate the sensitivity of the 
measured fluxes to the IMF parameters and present a bench- 
mark stellar population model for interpreting the results. 
Section [4] shows our primary result - that there is a strong 
correlation of the Ho to UV flux ratio with surface brightness 
of both Ha and R band emission, as well as weaker correla- 
tions of this ratio with other quantities. In Sec. [5] we discuss 
other physical parameters that can affect this ratio and argue 
that the IMF parameters are most likely driving the observed 
correlations. Section |6]places our results in the context other 
work, presents a physical scenario for pressure driven IMF 
variations driving the observed correlations, and discusses the 
implications of our results. Section [7] summarizes our results 
and suggests future observations. 

2. OBSERVATIONS AND MEASUREMENTS 

2.1. Sample Selection 

Our sample consists of 103 galaxies with both Ha obser- 
vations from SINGG and UV observations from SUNGG. 
The ultimate parent of b oth these surveys is the H I Parkes 
All Sk y Survey, HiPASS dMever et al.ll2004tlKoribalski et al.1 
120041) . M06 discuss the selection of the SINGG sample in 
detail. In brief, the sample was selected uniformly in the log 
of H I mass, Aim, with the nearest galaxies preferred in or- 
der to yield the best spatial resolution. Only H I properties 
were used to select the SINGG sample. SUNGG is a Cy- 
cle 1 GALEX Legacy Su rvey, whose targets were all c hosen 
from the SINGG sample (Wong 2007t lWong. et al.l2008l) . Be- 
cause of the large GALEX field of view (1.1 "diameter), some 
SINGG galaxies not in the original SUNGG selection are also 
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FIG. 1 . — The adopted Ha and FUV dust absorption in magnitudes, An a 
and Apuv, respectively plotted against each other. These are determined from 
the R band absolute magnitude (M' R ) and the FUV to NUV color (FUV - 
NUV)' corrected for Galactic dust absorption but not internal dust absorption. 
These quantities are shown on the alternate axes (top and right, respectively). 
The thin solid line shows the unity relation (Apuv —An a ). The heavy solid 
line sh ows the average ratio for our sample. The Apu y /Ahq ratio expected 
for the ICalzetti. et alJ2000l) and (Cardelli et al. 1989) dust absoprption laws 
are shown with the long dashed and dash - triple dot lines, respectively; while 
the the ratio Apu v /(2Ahci ) is shown by the short dashed and dash-dot lines 
respectively. 

included in the measurements presented here. As shown by 
M06, HiPASS sources often correspond to multiple emis- 
sion line galaxies (ELGs) in our narrow-band images (M06). 
By default we consider them part of the HiPASS source tar- 
geted. However, we have excluded five sources in multiple 
ELG fields from the analysis after inspecting the optical and 
UV images. Three of these are barely larger than the optical 
PSF and may be background galaxies (e.g. [Olll] emitters at 
z = 0.3), one may be part of a larger galaxy and not a galaxy in 
its own right, and the last is not totally within the CCD frame. 

2.2. Optical Observations 

The Ha and R band observations were obtained primarily 
with the CTIO 1.5m telescope, while some targets were ob- 
served with the CTIO 0.9m telescope to a similar depth, and 
at a similar plate scale. M06 describes the method of data ac- 
quisition, reduction and analysis. The data we use here are 
measurements of integrated flux and surface brightness. As 
detailed in M06 these are measured using concentric ellipti- 
cal aperture photometry after masking out obvious foreground 
stars, background galaxies and other image blemishes from 
the sky subtracted images. Plots of enclosed flux as a func- 
tion of aperture major axis radius are used to determine the 
total flux, F; the radii enclosing 50% and 90% of the flux, 7-50 
and rgo; and the effective surface brightness £ which is the 
face-on surface brightness within r^y. £ = F/(27rr| ), where 
here we use F to denote fluxes, or count rates in general. Later 
we become more specific and use F to denote integrated line 
fluxes and / to distinguish continuum flux density. 

The Ha signal to noise ratio is low, S/N < 2 within the out- 
ermost measurement aperture for 19 of the ELGs. In most of 
these Ha is clearly detected, but its signal is washed out over 
the large apertures needed to measure the total flux. For these 
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cases we estimate the total Fu a as twice the Ha flux within rso 
determined from the R band. This correction factor is consis- 
tent with the observation that the Ha and R emission typically 
have equal rso values (H06). Since the measurement errors 
are much smaller within rso this provides significant detec- 
tions for 5 of the low S/N sources. We use upper limits based 
on the r^o extrapolated flux for the remaining 14 sources. Of 
these, only one, HIPASS J 132 1-31, has no Ha emission dis- 
cernible in our images. The optical counterpart is a Low Sur- 
face Brightness (LSB) galaxy identified using H I synthesis 
imaging observations presented by R yan- Weber e t al. ( 2003; 
also shown by Grossi et al. 2007). IPritzl et all (120031) also 
note that no Ha is seen in their WIYN 3.5m observations of 
J 132 1-31. This source is discussed further in Sec 15.41 

2.3. Ultraviolet Observations 

The FUV and NUV GALEX images were typically 1500s 
(one eclipse) in dura tion, the same depth as the GALEX 
Nearby Galaxy Atlas dGil de Paz et alj|2007h . We make mea- 
surements from the combined flux calib rated images of each 
field taken from the GALEX pipeline (iMorrissev et al.ll2007l) . 
We perform our own background subtraction since the angular 
diameter of our galaxies is often > 3.2' whic h is the sky grid 
size u sed for GALEX pipeline processing (Mo rrissev et al.l 
120071) . Fluxes, radii (rso, rgo), and surface brightness in the 
FUV and NUV are measured in a manner analogous to that 
done for SINGG - using enclosed flux curves after masking 
out foreground and background objects. 

The apertures used by SINGG and SUNGG are configured 
independently. Hence the UV and optical apertures gener- 
ally are not matched, with the UV aperture almost always 
being larger. Comparison of the aperture areas with the rso 
and measurements, shows that the smaller of the UV or 
optical aperture should recover the majority of the flux at all 
wavelengths in almost all cases. Specifically, if the Ha and 
UV light trace each other (as expected for a universal IMF) 
the bias induced by mismatched apertures in \og(Fn a / /fw) 
should be < 0.05 dex in 81% of the sample, may be as high 
as 0.3 dex in 17% and may be higher in 2% of the sample 
(2 cases). The actual bias depends on how much light is out- 
side the optical measurement aperture, and may be consid- 
erably lower than these estimates if the IMF at large radii is 
deficient in high mass stars as found in some spiral galaxies 
dThilker et al.ll200HlBoissier et al.ll2007l) . 

2.4. Error bars 

The error bars we display account for the measurement er- 
rors only. The errors in the optical quantities are discussed 
in detail in M06, and include terms for the sky uncertainty 
and the continuum subtraction. The errors in the UV quan- 
tities includ e terms for photon statistics an d the background 
uncertainty dWongll2007i IWong. et alj|2008l) . 

2.5. Corrections for Emission Line contamination and 
Galactic Dust Absorption 

A variety of corrections to the measurements are required to 
transform them to intrinsic quantities. M06 details the correc- 
tion of the Ha measurements for [N II] emission and Balmer 
absorption line contamination. The Ha, R, FUV, NUV mea- 
surements are all corrected for foregro und Galactic dust ab- 
sorption using the redde ning maps of ISchlegel et al.l (1 19981) 
and the extinction law of Card elli et al.l u 198S 
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FIG. 2. — Sequences of population synthesis models showing the effects 
on the luminosity in Ha and FUV of varying M u at fixed 7, relative to 
a fiducial model having M u = IOOMq. Three sequnces corresponding to 
7 = -3.3,-2.35,-1.3 (differentiated by line style as shown in the legend) 
are shown for each band with the Ha sequences shown as (red) connected 
dots and the FUV sequences as (purple) connected triangles. Horizontal lines 
from top to bottom indicate luminosities of 95%, 50%, and 5% relative to the 
fiducial model. See the text for further details of the models. 

2.6. Correction for Internal Dust Absorption 

Internal dust absorption can have a substantial effect on in- 
cident fluxes, especially in the UV, causing systematic biases 
in ^ho/Zfuv- We correct the optical and UV data using sepa- 
rate empirical correlations. 

For the UV data we base the correction on a fit to the corre- 
lation between the ratio of total far infrared (FIR) flux and 
the FUV flux and the FUV - NUV color. The total FIR 
flux w as determined from IRAS fluxes following iDale et afl 
d2001l) for 73 galaxies in the SUNGG sample that were de- 
tected by IR AS. The IRX-/3 d ata and fit for our sample 
is shown in IWong. et al.l (120081). Similar cor r elation s are 
shown by ICortese et alTd2006l) . iBoissier et ail d2007l) , and 
iGil de Paz et al.l d2007l) . The fitted relationship thus provides 
a method to recover the light absorbed by dust and re-radiated 
in the far-infrared. 

The optical fluxes were corrected for internal dust us- 
ing a relationship between Ha dust attenuation derived 
from Balmer line ratios and R band absolute magnitude 
prior to i nternal dust cor re ction, M' R , found and pre- 
sented by Helmbo ldt et al.l d2004l) us i ng da t a from the 
Nearb y Field Galaxy Survey lJansenl l2000t iJansen et al.1 
2000). The R band internal dust attenuation, in mag- 
nitudes was taken to be half that of Ha to account for 
"differential" dust absorption: emission lines are more ef- 
fected by dust than the stellar continuum (Panelli et aD 19881: 
ICalzetti. Kinney & Storchi-Bergmanlll994l) . presumably be- 
cause high-mass star formation is more highly correlated with 
dust than the general field star population. M06 present a pop- 
ulation synthesis model attenuated by dust using this prescrip- 
tion in order to explain the relationship between Lr and the 
ratio of Ha to FIR fluxes, fka/^FiR- The model used had 
a solar metallicity, a Salpeter IMF over the mass range of 
0.1 to 100 Mq, and a constant SFR f or a duration 100 Myr , 
and was calculated using Starburst99 dLeitherer et al.| [l999). 
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TABLE 1 

Contribution to luminosity by stellar mass 



contribution 




Ha 






FUV 




level 


7 = -3.3 


7 = -2.35 


7 = -1.3 


7 = -3.3 


7 = -2.35 


7 = -1.3 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


5% 


17 


23 


28 


2.5 


3.3 


7.7 


50% 


40 


54 


71 


5.7 


16 


47 


95% 


91 


95 


97 


35 


78 


94 



NOTE. — The values in this table correspond to the intersection of the horizontal lines 
with the curves shown in Fig. [2] They show the contribution of different stellar masses 
to the luminosity of a solar metallicity stellar population that has been forming stars for 
1 Gyr. The stellar population has a single power law IMF with mass range of 0.1 to 100 
Mq. The table gives the upper end of the mass range, in Mq, that contributes the first 
N% (specified in column 1) to the luminosity in Ha (columns 2-4) and FUV (columns 5 
- 7) for three different IMF slopes 7. Thus for a Salpeter IMF (7 = -2.35) 5%, 50%, and 
95% of the Ha luminosity comes from stars having masses up to 23 Mq, 54 A4q, and 
95 Mq respectively, while 5%, 50%, and 95% of the FUV luminosity comes from stars 
having masses up to 3.3 A4q, 16 Mq, and 78 Mq respectively. 



while the dust followed the Calzet tiTet al.l (|2000) dust attenu- 
ation model. While this model accounted for the shape of the 
Fna/Fpm versus L« relationship, there is an offset: Fn a /^fir 
is predicted to be too high by a factor of ~ 3. There are vari- 
ous possible explanations for this discrepancy, some of which 
are discussed by M06. In particular, they note that adopting a 
steeper 7 and/or lower M. u will help alleviate the discrepancy. 
Hence, the low measured Fjj a /Fp[R values may be related to 
the low Fh q //fuv values that are the primary concern of this 
investigation. Other possible factors that could contribute to 
the M06 Fhq/Ffir model being too high include: too short of 
a duration for star formation, a metallicity that is too low, es- 
caping ionizing photons, choice of stellar population models, 
a dust absorption model that is too differential, and an error 
in the FIR bolometric correction. It is beyond the scope of 
this paper to improve the Fhq/Ffir versus M' R model of M06. 
However, except for the dust bolometric correction, we do 
address all of these issues here with regard to our Fn a /fpw 
results. Since the overall shape of the M' R - Fn a J -Ffir relation- 
ship agrees with the M06 model, then using this optical cor- 
rection should yield the correct Ha and R band fluxes modulo 
a zeropoint offset. 

Figure[T]compares the internal dust absorption in the FUV, 
Afuv, with that in Ha, An a , derived from the relations de- 
scribed above. It demonstrates that the internal dust correc- 
tions are broadly consistent with vectors from observed red- 
dening laws, especially those that account for the differential 
extinc tion between emiss ion lines and the stellar continuum 
(e.g. ICalzetti7etaill2000l) . Since these are transformations 
from the UV color and M' R respectively, these quantities are 
shown on the alternate y and x axes. For our adopted dust 
absorption law, the x axis also shows 2Ar where Ar is the R 
band absorption. The Pearson correlation coefficient r xy be- 
tween Apuv and An a is a 0.59. The probability to achieve this 
at random from uncorrelated quantities is 4 x 10~ 12 . The aver- 
age ratio, (Apuv /Au a ) = 1 -68, is plotted as a thick solid line in 
Fig-D The da shed lines show that this is close to what is ex- 
pected for the Calze tti. et al.l (|2000) starburst dust absorption 
model, which includes differential absorption, applied to pop- 
ulation synthesis models (as described below). This model 
results in A FU v/A Hct = 1.54 and A FUV /(2A R ) = 1.50. The ob- 
served average ratio is somewhat shallow for Ha compared to 
"normal" dust reddening laws, where dust and gas are treated 
the same. For example, usin g the same modeling we calcu- 
late A FUV /A Hct = 2.50 for the lCardelli eTail (119891) reddening 



law (shown as the dash - triple dot line in Fig. [TJ. A better 
agr eement is found with the continuum, Afuv/(2A«) = 1.60 
for lCardelli etail d!989t) (dash - dot line). 

The amplitude of the dust corrections are also consistent 
with what is known about normal star forming galaxies. Our 
sample has a range of Au a from 0.13 to 1.8 mag with an 
average and median of 0.71 and 0.56 mag respectively. We 
compare this to the Spitzer Infrared Nearby Galaxy Survey 
(SINGS) for which galaxy averaged An a can be estimated 
from a weighted average of Ha and Spitzer 24/xm fluxes , 
following the calibra tion presented by Cal zetti et al.l (120071) . 
iKennicutt et al.l d2008l) do this and find the SINGS sample has 
Ana ranging from to 2.5 mag with the median at ~ 0.7 mag. 
That is, broadly consistent with what we find. For the FUV, 
we find Apuv ranging from to 3.03 mag, with an an aver- 
age an d median of 1.2 mag. From Fig. 8 of Gil de P az et al.l 
(2007) it can be seen that Apuv typically lies between 0.6 and 
2 mag in their angular diameter selected sample. Despite the 
general consistency of our dust corrections with expectations, 
one must always be careful of the biases dust may cause, es- 
pecially in cases like ours where the estimate of the correc- 
tion is indirect. For example, Fig.[T]shows that relative to the 
(Afuv /An a ) line, galaxies with low An a (< 0.6 mag) tend to 
have high Afuv, while for high An a (^ 1 mag) Afuv tends 
to be low compared to An a - In Sec. 15.11 we show that this 
difference does not have a significant effect on our results. 

2.7. Quantities and units 

The primary quantity we work with here is the ratio of 
the Ha line flux Fn a , to the FUV flux density /fuv (inci- 
dent flux per unit wavelength), Fn a /fFuv, which has units 
of A. We compare this to the effective surface brightness 
in Ha and the R band, £na and £#, respectively. We 
show the Ha surface brightness Sea in units of W Kpc~ 2 (= 
5.1 x 10 47 ergcm~ 2 s _1 arcsec" 2 ) and the R band surface bright- 
ness in units of Lr q Kpc" 2 , where the sun's R luminosity is 
L r .q = 4.39 x 10 22 W A" 1 {M r .q = 4.61 ABmag). 

3. STELLAR POPULATION MODELS 

3.1. Sensitivity to IMF parameters 

The luminosity in Ha and FUV of a stellar population 
arise from different but overlapping mass ranges of stars; 
this drives the sensitivity of Fhq.//fuv to the properties of 
the upper end of the IMF. This is illustrated in Fig. [2] which 
shows the effect of varying Ai„ on the expected Ha and 
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FIG. 3. — The ratio of Ha line flux to FUV flux-density as a function of (a) the face-on Ha effective surface brightness, and (b) (right) the face-on R band 
effective surface brightness £/;. All quantities have been corrected for Galactic and internal dust absorption as described in the text. The solid line shows the 
iteratively clipped least squares fit to the data, while the dashed lines show the final clipping limits. The dot-dash line shows the Fhq//fuv ex pected for our 
fiducial stellar population (eq.f5J. The vectors i n the lower right porti on of the panels shows the effect of internal dust absorption according to the Calz etti. et alj 
12000) starburst attenuation law (red vector) the Cardelli et al. 1 1989) Milky Way dust extinction law (orange vector), as well as our method (black vector). The 
length of the vectors in log(/"Ha//FUv) in all cases is set to the average correction we deduce. Galaxies corresponding to measurements A and B are shown in 

Fig.m 

luminosity is more sensitive to 7 than M. u . 

3.2. Fiducial stellar population model 

From the models described above, we adopt the model with 
7 = -2.35 and Al,, = 100 Mq as our fiducial stellar population 
model. These parameters, or similar values are often used 
to chara cterize normal sta r forming stellar populations. For 
example iKennicutl (1 19981) adopted these parameters for his 
SFR calibrations, alt hough they are base d on different stellar 
population models (Mad au. et al.H l998). From our fiducial 
model we d erive the followin g SFR calibrations, and compare 
them to the iKennicuttl (TT998 ) calibrations (in parenthesis): 



FUV luminosities for stellar populations having a single slope 
IMF with Mi = 0.1 M®, and a constant fixed SFR for the 
stars having M.+ < The stellar population synthesis 

models w ere calculated using t he online Starburst99 pack - 
age (v5.1. lLeitherer et alj[l999l: IVazquez & Leifhererll2005l) . 
They use solar metal li city Padova g r oup evo lut ionary tracks 
dBressan et all Il993t iFagotto et all | 1994allbl: iGirardi et alj 
and model at mospheres from Pauldr ach et al.l |2001) 
and iHillier & Millen (119981) : these are the default options in 
Starburst99. Case B recombination was assumed to con- 
vert the output o f ionizing photons to an Ha luminosity 
dOsterbrockll 19891) . The FUV luminosity was determined by 
applying the GALEX FUV transmission curve to the spectra 
generated by Starburst99. We adopt a fixed duration At = 1 
Gyr for the star formation so the births of O and B stars are 
in equilibrium with the deaths. Figure |2] shows that only stars 
with Al* > 3A4® contribute significantly to the Ha or FUV 
lumin osities and these have main sequence lifetim es < 325 
Mvr (ISchaerer et aljfl993afbl: IFagotto et alj|1994albb . Three 
model sequences with IMF slope 7 = -1.3,-2.35, -3.3 are 
illustra ted. The central value corresponds to the [Salpeter 
(1 19551) slope which is consistent with more recent determi- 
nations of Kroupa (2001, 2002). while the other two val- 
ues are meant to represent shallow and steep extremes to 
the slope. The y axis plots the luminosities relative to the 
case were M. u = lOCM©. This corresponds to the highest 
stellar mass determined dynamically us i ng binary star orbits 
dBonanos etal.ll2004t iRauw et al.ll2004t ISchnurr et alll2008l) 
and is a typical A4„ adopted in the literature. The horizontal 
lines indicate the masses below which the contribution to the 
luminosities is 5%, 50% and 95% of the total luminosity if the 
true M. u is 100 M.®. We tabulate the values correspond- 
ing to these contributions in Table Q] 

The Ha curves are very steep and closely spaced together. 
This indicates that the mass range contributing to the Ha lu- 
minosity is more sensitive to M.„ than 7. The opposite is true 
for the FUV luminosity - the curves are relatively shallow for 
near M. u and the shapes are very different as a function 
of 7 this indicates that the mass range contributing to the FUV 



SFR(Ha) _ L Ha 
IMqvt 1 ~ 1.04(1.27) x 10 34 W' 



SFR(FUV) 
lAloyr" 1 



9.12(9.09) x 10 32 W" 



(1) 



(2) 



Here the UV calibration of iKennicuttl (fl998h at A = 2150A 
was transformed to the pivot wavelength Xp = 1535A of the 
GALEX FUV filter assuming the an intrinsic power law FUV 
spectrum fx oc A' 3 with UV spectr al slope J3 = -2 . To com- 
pare results to those that adopt the Kroupa (2001) IMF (e.g. 
iBrinchmann et alJl2004fclKauffmann et al.ll2003h. the SFR es- 
timat es here should be divided by 1 .5 dBrinchma nn et all 
l2004h . 

The ratio of eq.|2]and eq.Q] gives Fh q //fuv f° r me fiducial 
model: 



/fov 



= 11.3 



(3) 



We show this ratio in various plots as an indication of the a 
priori expected Fna/fpuv- 

In order to test the sensitivity of our results to the choice of 
stellar population mo dels we also calculated models us i ng the 
codes PEGASE ( v2: iFioc & Rocca-yolmerange!ll997ll 19991) 
and GALAXEV dBruzual & Charlotll2003h . In both cases we 
adopted an IMF, metallicity, and SFH identical to the fidu- 
cial stellar population. Both the PEGASE and GALAXEV 
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TABLE 2 
Properties of example galaxies 



Parameter 


Example A 


Example B 


units 


Name 


UGCA044 


NGC1566 




HIPASS+ 


J0249-02 


J0419-54 




Morphology 


IB(s)m: 


(R',)SAB(rs)bc 




logCMm) 


8.85 


10.19 


Mq 


log(Z*) 


8.29 ±0.04 


11. 09 ±0.01 


Lr,q 


M FUV 


-14.76 ±0.04 


-20.66 ±0.00 


ABmag 


rso 


2.37 ±0.05 


6.76 ±0.09 


Kpc 


log(S Ha ) 


30.82 ±0.07 


32.89 ±0.05 


WKpc" 2 


log(£ R ) 


6.87 ±0.01 


8.70±0.01 


Lr.q Kpc' 2 


log(F Ha //Fuv) 


0.71 ±0.08 


1.40 ±0.05 


A 



FIG. 4. — Example galaxy A - UGCA44 (HlPASS J0249— 02) shown in Ho, R band, and FUV in red, green, and blue respectively in the left panel, while the 
right panel shows the R band only in an inverted gray scale. The three images that comprise the right panel have been convolved to have the GALEX resolution 
of 6.5", while the image on the left panel is shown at the original resolution of the SINGG data (typically ~ 1.6"). The scalebar in the right panel has a length 
of 10 Kpc. 

Figure [3] shows our primary result: a strong correlation be- 
tween Fua/fFW and surface brightness in both Ha (panel a), 
and in the R band, £r (panel b). One striking aspect of the 
Figure is that more than half the sample has ■Fhc<//fuv below 
normal. Except for possibly one weak Ha detection, and one 
galaxy with an upper limit to Fu a / /fu v , this includes all of the 
LSB galaxies having log(S HQ ) < 31.7 or log(S s ) < 7.5. Or- 
dinary least squares bisector fits were performed to the data, 
with iterative clipping of the points having residuals outside 
of ±2.5<7 V from the best fit line yielding 

log(f^) = (-12.75 ± 1.02) + (0.43 ±0.03)log(S]H Q ), (4) 

log(f^) = (-3.16 ± 0.28) + (0.52 ± 0.03) log(E*). (5) 

These are shown as solid lines in Fig.[3]while the clipping lim- 
its are shown as dashed lines. Data with only Ha upper limits 
were not included in any of our fits. The quality of the corre- 
lations and fits are very similar with each having a Pearson's 
correlation coefficient r xy = 0.65, a probability of ~ 4 x 10~ 12 
that the correlation could have occurred from random data, 
and an rms dispersion of 0.22 dex about the log(fHo,//Fuv) 
residuals (after clipping). 

Optical surface brightnesses are not the only quantities that 
correlate with Fhq,//fuv- This is illustrated in two ways. First, 
in Fig.|4]we show example images of galaxies at either end of 
the correlations as indicated in Fig.[3j some properties of these 
two examples are given in Table[2] This illustrates the correla- 
tion with morphology along this sequence. The LSB example 
galaxy, UGCA44 is a dwarf irregular galaxy without coher- 
ent structure. The High Surface Brightness (HSB) example, 
NGC 1566, the brightest member of the Dorado group, is a 
strong regular two arm spiral galaxy. These two examples 
are typical of the morphological variation along the sequence, 
and consistent with what is known about the sequence of star 
forming galaxies: LSB sources are typically low luminosity 
irregular systems, HSB galaxies tend to be high luminosity 
spirals (as well as some blue compact dwarfs), galaxies with 
moderate surface brightness are a mixture of late type spirals 
and brighter irregulars. 

Second, Fua/fFuv correlates with other global quantities 
such as total luminosity L K and rotational velocity V rot (deter- 
mined from the H I line widths as outlined in H06), which is 



models use the same Padova group sources for their evolu- 
tionary tracks as does our Starburst99 model. The differ- 
ences are in the Stellar atmosph eres. The PEGASE model 
uses a library of observed stars dFioc & Rocca-V olmerange 
1 1997b while GALAXEV us es various BaSeL t h eoreti cal at- 
mospheres as described by iBruzual & Chariot d2003l) . We 
used our own software to calculate / FUV as described above, 
while the models provide either L# Q (PEGASE) or ionizing 
photon flux from which we calculate L# Q (GALAXEV). We 
calculate an equilibrium /*hq//fuv = 10.8, 11.8 for the PE- 
GASE and GALAXEV models respectively. Hence, differ- 
ences in stellar population models can affect the Fhq//fuv 
calculations at the ~10% level. Metallicity also effects the 
model results; metallicity sensitivity is considered in Sec. 15.51 
The SFR conversion factors and equilibrium Fhq//fuv val- 
ues all assume that the IMF is universal. However, as shown 
here, it is likely that the IMF is v ariable and often is not 
consis tent with our fiducial model. Pflamm-Altenburg et al. 
(2007) examine the implications to Ha based SFR estimates 
from an IGIMF that depends on the SFR, as one would expect 
for star formation dominated by clusters of finite mass. Then 
the conversion factors can vary from the above by orders of 
magnitude. This scenario is discussed further in Sec. 15. 41 and 
SecIO 

4. RESULTS 
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FIG. 4. — continued. Example galaxy B - NGC1566 (HlPASS J0419— 54). The right and left panels have the same display levels as the other example. The 
projected area of these images is larger, while the scalebar length remains 10 Kpc. 



shown in Fig. [5] Only HlPASS sources corresponding to sin- 
gle ELGs with optical axis ratio a/b> 1.4 are plotted in panel 
b, in order to limit the plot to galaxies likely to be dominated 
by rotation. Using the same fitting algorithm as described 
above, we find 

log(^) = (-1 .64 ±0.15) + (0.27 ± 0.01)log(L s ), (6) 
log(^) = (-0.97 ±0.17) + (0.96 ±0.09)log(V rot ). (7) 

These correlations are weaker than those shown in Fig. [3] hav- 
ing correlation coefficients of 0.61 (0.57) while the (clipped) 
rms scatter in the log(Fn a / /fuv) residuals is 0.19 (0.32) dex 
in panel a (b). Nevertheless, the probability that these cor- 
relations could have come from randomly drawn data is still 
low: 10~ 10 for the correlation with log(L«), and 3 x 10~ 4 for 
the correlation with log(V rot ) (low, but much higher than the 
other correlations shown here because of the smaller sample 
of elongated single ELGs). 

The range of luminosity and mass in these correlations are 
too large for galaxies to evolve from one end of the correlation 
to the other. Instead, Fh q //fuv is one of the global properties 
(although not a fundamental one) that vary regularly along 
the sequence of late type galaxies. O ther well know n correla- 
tions along this sequence include the Hubble; (1 19261) morpho- 
logical sequence, th e Tully-Fisher (1977) relation, the mass- 
metallicity relation (Tremonti et al. 2004| ) and the Universa l 
Rotation Curve dPersic & Saluccil 1 1 99 It iPersic et ai1ll996l) . 
The main physical parameter that drives most of these cor- 
relation appears to be halo-mass. However, the correlations in 
Fig.|3]are stronger than those in Fig.|5]indicating that star for- 
mation intensity and stellar mass density are more important 
for driving the Fh q //fuv correlations than halo properties. 

Both the range of /*hq,//fuv values, covering over a factor 
of ten, and the correlation with surface brightness have major 
implications, as discussed in Sec. 16.31 

5. PARAMETERS AFFECTING F Ha //Fuv 

Physical properties that can affect Fho/Zfuv include dust 
absorption, the detailed SFH, the porosity of the interstel- 
lar medium, stochastic effects from the number of stars, the 
metallicity of the stellar populations, and the properties of the 
IMF. We now consider, in turn, which of these could drive the 
observed correlations, saving the IMF for last because it is the 
only option we can not rule out. 



5.1. Dust 

The dust corrections are too small and in the wrong sense 
to have spuriously caused the observed correlations. More- 
over, even without any dust corrections, a large fraction of the 
sample have too few ionizations for the observed UV light 
compared to expectations for normal star forming popula- 
tions. These points are illustrated in Fig. [6] which is simi- 
lar to Fig. [3] but here the quantities are shown before inter- 
nal dust correction (and denoted with a prime: E' Ha , and 
^Ho/Zfuv)- We see strong correlations between the quanti- 
ties even before dust correction. Here r xy = 0.66 (0.69) for the 
correlation between /^ Q // FUV and 2^(1]^), while the proba- 
bility that these are spurious correlations from random data is 
7 x 10" 13 (3 x 10" 14 ). The best fits to these data are 

log (%*■) =(-11.32±l. 14) + (0.39 ± 0.04) logCE^J, (8) 

\ J FUV / 

logfllM =(-3.36±0.32) + (0.58±0.04)log(E'). (9) 

\ J FUV / 

As before, the fits are OLS bisector fits with an iterative 
2.5cr clipping. The scatter about the fit is 0.22 (0.23) dex in 
log(/^ a //p UV ) for the E[j Q (E^) data, nearly identical to the 
scatter seen in Fig. [3] 

There are several reasons why reddening correction does 
not greatly effect the results. First, the average dust correc- 
tion, shown by the black arrows in Figs. [3] and [6] is not large 
compared to the range of the correlations. Second, The differ- 
ence in the dust correction between the high and low surface 
brightness ends is also not that large. This is shown by di- 
viding the sample into quartiles of F^ a /fp UV . The quartile 
boundaries are at log(F^ a / / FUV ) = 0.91, 1.04, and 1.26 and 
shown with dotted lines in Fig. [6] The average of the uncor- 
rected points within the quartiles are shown as large pink tri- 
angles, while the average of the same data points after dust 
correction is shown with large cyan diamonds; thick lines 
show OLS-bisector fits to these average points. Systematic 
changes in the Apuv /Ah q ratio noted in 32.6l are also not large 
enough to seriously effect the results. 

Figure [6] also shows that nearly half the sample have 
Fn a / /fuv below the value for our fiducial stellar population 
model. The first quartile has <log(/^ Q // FUV )) = 0.73 ±0.10, 
where here the error is the average error on log(/^ Q // FUV ). 
Hence, the galaxies in this quartile have log(/^ Q // FUV ) val- 
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ues that are on average 3a lower than the fiducial model even 
before dust correction. 

Could it be that we are under correcting for dust absorp- 
tion? This seems unlikely, because the amount of dust needed 
to flatten the relationships shown in Fig. [3] is implausibly high 
and would only exacerbate the problem of low ^hq//fuv val- 
ues. To flatten the relationship we consider the amount of dust 
required to bring the quartile (log(f^ Q ,//p UV )) values equal to 
that of the first quartile. In other words, this assumes that 
the first quartile is virtually unaffected by dust. Since dust 
almost certainly effects all measurements this will produce 
an underestimate of the dust correction needed. In the sec- 
ond, third, and fourth quartiles (log^a/ZFUv)) = °- 97 ' !- 15 ' 
1 .43 respectively. Using (A F uy/A Ha } = 1.68 (Sec. 12.6b then to 
lower these to the value observed in the first quartile requires 
Ana = 0.9 1.6, and 2.6 mag, respectively, and Afuv = 1-5, 2.6, 
and 4.4 mag, respectively. This is equivalent to an overall av- 
erage Ahq = 1.3 mag and Afuv = 2.1 mag, or about twice the 
value estimated in Sec. 12.61 

Much of the scatter in the observed correlations may come 
from uncertainties in the internal dust correc ti on. From 
the scatter in the IRX-/3 fits of IWong, et all d2008l) and 
ICortese et al.l (12006) combined with th e relationsh i p betw een 
IRX and FUV dust absorption A F uv of iBuat et all (120051) we 
determine that about 0. 14 dex uncertainty in /fuv results from 
the IRX-/3 relation. The correction to Fa a for dust absorp- 
tion is tied to Balmer Decrement measurements. iKewley et al.1 
(2002) show that after dust correction the SFR estimated from 
Ha and FIR emission agree well in the calibrating Nearby 
Field Galaxy Sample. Using the online version of their table 
1, we find a dispersion in their mean log(SFRpiR/SFRHa) of 
0.16 dex (after an iterative 2.5a clipping). Assuming this is 
the intrinsic scatter due to dust correction affecting Fn a , and 
that it is uncorrected with the scatter in /fuv, then the uncer- 
tainty in log(,pHa//Fuv) due to dust may be as high as 0.21 
dex. This is very close to the observed scatter in Figures [3] 
and [5] However, this might be an overestimate of the scatter 
due to dust since SFRh q estimates of Kewley et al. are also 
affected by the IMF. 

Driver et al. (2007,2008) find strong inclination dependent 
differences in the luminosity functions of the bulge and disk 
components of galaxies which they attribute to dust absorp- 
tion. While inclination effects may be seen in, and removed 



from, the UV with our color based Afuv estimate, our Ah q 
estimate is just based on luminosity and any inclination term 
should remain and effect or Fh q //fuv estimates. However, 
we found insignificant correlation coefficients r xy = 0.13, 0.01 
for the residuals of the fits in Fig.[3]panel (a) and (b), respec- 
tively, and the axial ratio (a/b) of the sources in the optical, 
demonstrating that inclination induced dust absorption is not 
significantly contributing to the scatter. Further work is re- 
quired to determine what is causing the scatter in these rela- 
tionships. 

5.2. Star Formation History 

Sharp changes in the SFR can strongly effect Fh q //fuv due 
to the different lifetimes of O and B stars. To explore the ef- 
fects of a variable SFR on our correlation we built photomet- 
ric models of a temporary increase, or "burst", and decrease, 
or "gasp", on an otherwise constant SFR population buildup. 
The models were created using Starburst 99 to calculate the 
evolution of the spectrum of a simple stellar population (i.e. 
instantaneous starburst). The evolution of the relevant lumi- 
nosities (Ha, FUV and R) for an arbitrary SFR as a function 
of time, t , were calculated from the simple stellar population 
results using a separate program. These models employed the 
same fiducial IMF outlined in Sec. 13.21 

Figure [7] shows the burst and gasp SFHs that we modeled, 
the evolution of the fluxes Fn a , /fuv, and /r, and the evolu- 
tion of the ratio -Fho//fuv- The models are relative in the 
sense that the absolute SFR is not important and is set by 
adding appropriate zeropoints. Hence the SFR and fluxes are 
shown relative to those before the burst or gasp. The func- 
tional form that we adopt is a constant SFR with a Gaussian 
enhancement or depression representing a burst and gasp, re- 
spectively. The model parameters are the FWHM duration, 
At, of the event, the ratio of maximum to minimum SFR, A, 
and the time to of the center of the event. All models have a 
finite base SFR. We do not allow the SFR to completely turn- 
off because our sample shows that almost all late type galaxies 
have at least some recent high mass star formation (M06). For 
the burst models, the maximum SFR occurs during the event 
and the minimum is the base level, while for gasp models, the 
minimum SFR occurs during the event and the maximum is 
the base level. All models have to = 10 Gyr. Absolute time is 
important when considering R band fluxes; our choice places 
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FIG. 6. — The quantities shown in Fig.fJ]are plotted here prior to dust absorption correction, and hence are denoted with prime (') symbols. The dotted lines 
separate the sample into quartiles of log(/^ Q //p UV ). The pink triangles show the average qauntities of the data points in each quartile, while the blue diamonds 
shows the average of the same points after dust correction. The thick pink and cyan lines are least squares fits to the respective average points. The meaning of 
the other symbols and lines are the same as in Fig. [5] 



the burst or gasp at roughly the Hubble time. 

We explored a range of model parameters and found that the 
salient features of these models can be illustrated with models 
having At = 10, 100, 1000 Myr, and A = 2, 10, 100. Figure 
shows that very large excursions in Fua/fpuw are possible, 
especially for large A and short Af . The Af = 1 Gyr models 
produce only very weak excursions in Fn a /fFuv because the 
timescale is long compared to the lifetimes of both the O and 
B stars. 

Figure [8] plots the tracks the models make in the 
Fna/fpuv versus Sh q and Sr planes. All models con- 
verge at our adopted surface-brightness zeropoints and 
log(FH Q //Fuv[]) = 1-055, the equilibrium value (eq. |3}. We 
chose log(E Ha [WKpc- 2 ]) = 32.3 and log(E^[L RiQ Kpc- 2 ]) = 
8.2 as the zeropoints corresponding to the surface brightness 
immediately prior to the burst or gasp (f = fo-2Af ) since these 
values are near the middle of the pack with respect to other 
galaxies having similar Fn a / /fuv- We also ran the models 
with the same IMF, metallicity and SFHs but based on PE- 
GASE model calculations. Other than the zero-point offset 
noted in Sec. 13.21 these models make indistinguishable tracks 
from those shown in Fig. [8] 

Both the burst and gasp models can produce large corre- 
lated excursions in Eh q and Fhq//fuv mimicking the ob- 
served correlations. This is especially true for the gasp mod- 
els. However, in cases of large A and short Af the tracks can 
enter unpopulated regions of the ^hq/Zfuv versus £hq dia- 
gram during the post-burst phase as Fjj a /fyiyv drops and Eh q 
fades to its original value. Fine tuning of the SFH to have a 
longer duration to the turn-off phase may alleviate this prob- 
lem. The Fft a /fp\jy versus plane is more difficult to ex- 
plain with burst or gasp models. This is because short Af 
events do not last long enough to significantly change 
while the long St events that can cause large excursions 
have little effect on Fn a /fpuy. Thus a SFH including a recent 
burst or gasp can not simultaneously account for the correla- 
tions of Fn a /fyuv with both and Eh«- If bursts and gasps 
cause the spread of -Fife //fuv then the SFH of galaxies must 
be synchronized so that LSB galaxies would all be in a gasp 
or post-burst phase. We can think of no physical mechanism 
that can naturally result in such a contrived scenario. 



While the burst and gasp models can produce large excur- 
sions in Fhq//fuv, this requires short Af and large A events. 
What amplitudes and durations are likely for intermittent star 
formation? Causality dictates that the event duration should 
not be shorter than the dynamical time, Td yn since it is the 
timescale for disturbances (such as the build-up or removal of 
ISM) to travel through a galaxy. Figure [9] displays the dis- 
tribution of T^yn calculated from the H I profile width for the 
galaxies in our sample with axial ratio a/b > 1.4. These cal- 
culations follow the method of H06 for calculating V rot H I 
line width and adopt Td yn = r or b /Vjtt, where r 01 -b = 2nr5o /Vrot- 
This definition effectively adopts a flat rotation curve which 
is evaluated at the Ha r^, and represents an average Td yn for 
high mass star formation in the galaxies. For rotation curves 
that are rising at rso, this definition underestimates Ta yn . From 
Fig IU it can be seen that the median Td yn ~ 60 Myr. Numeri- 
cal simulations of starbursts indicate durations as low as ~ 60 
Myr for the final sharp peak in the SFR during a merger, with 
the peak amplitudes on the o r der of 10-80 times the base level 
Mihos & Hernquisu(Tl994allbl[T996l) . These simulations show 
that enhanced star formation during a merger lasts several 100 
Myr. Mergers are dramatic rare events and our sample shows 
very few examples of these (< 10%, M06). Normal vari- 
ations in SFR should be weaker and longer. For example, 
morphological and SED analyses of starbursting dwarf galax- 
ies suggest A < 3 and Af in the range of se veral 100 Myr 
(iMeurer et al.ll 19921: IPapaderos et al.lll996aTlbt iMarlowe et all 
1999). Our models indicate that such events may add scat- 
ter to our correlations but can not explain the large range in 

■fka//FUV- 

5.3. Porosity of the Interstellar medium 

The observed Fn a /fpw correlations could also occur if the 
escape fraction of ionizing photons, / esc , inversely correlates 
with surface brightness. This explanation would require that 
the galaxies in the lowest quartile of /*h q //fuv are losing 
> 50% of their ionizing flux compared to our fiducial model. 
Escaping ionizing photons have been det ected in UV bright 
star forming galaxies at a redshift z ~ 3 (ISteidel et alJl200U 
Shapl ev et al.ll2006t llwata et al.ll2008l) . In the local universe, 
direct measurements of / e sc have been limited to HSB star- 
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FIG. 7. — Models of the photometric evolution of "burst" and "gasp" SFHs 
are shown in the left and right panel repectively. The top panels shows the 
relative star formation rate for the adopted Gaussian profile SFH. The x axis 
on all panels is the time relative to the event center to. Three FWHM dura- 
tions At = 10, 100, 1000 Myr are displayed using different colors: magenta, 
green, and black respectively for burst models; and red, cyan, and blue, re- 
spectively for gasp models. Three relative amplitudes, defined as SFR maxi- 
mum:minimum are shown - 2:1, 10:1, 100:1 and distinguished with increas- 
ing line thickness. The color and thickness coding are repeated in the re- 
maining panels and the next figure. The middle panel shows the evolution 
of luminosity relative to the pre-event luminosity (defined as tg - 2 At) in the 
Ha, FUV, and R band distinguished using solid, dashed and dotted lines re- 
spectively. The bottom panels show the evolution log(Fu a / /fuv)- 

burst galaxies and have mostly been upper limits in the range 
./• esc < 0.05-0. 1 (Leither er. et all 1 199 a iHurwitz et al J 119971: 
iDeharveng et al.l2001l) . The strongest c laim for a direc t detec- 
tion of /esc ~ 0.04-0.11 in Haro 1 1 bv lBergvall eTail d2006l) 
is disputed by iGrimes eTail d2007l) who find / esc < 0.02 for 
the same data. The best case for escaping ionizing photons 
from a normal galaxy comes from the very faint Ha emis- 
sion detected in the Magellanic Stream and High Velocity 
Clouds whi ch are thought to be partially ionized by the disk of 
our G alaxy (Bland-Hawthorn & Malone vll999t|Putman et ail 
2003). These measurements require about 6% of ionizing 
photons to escape normal to the Galactic plane, or, about 
/esc ~ 1% -2% when isotropitized over all angles. No di- 
rect measurements of / esc from LSB galaxies have been at- 
tempted. Hence we can not directly rule out the possibility 
that / esc inversely correlates with surface brightness. 

However, there are two strong arguments against this no- 
tion. First, it implies the existence of "naked" O stars which 
we should see directly in HST images of the nearest low sur- 
face brightness galaxies. Instead, the color magnitude dia- 
grams of dwarf irregulars are typically deficient in high mass 
stars ( e.g.lTosi et alj[l99lt iGreggio et aDll993t iMarconi et aD 
119951: lTolstovlll99aV While one may argue that the high 
mass stars are completely hidden by dust this would be in- 
consistent with them being naked. Furthermore, high mass 
stars can be isolated in large numbers from color-magnitude 
diagrams of higher surface-brightness dustier galaxies (e.g. 
Anger etti et ail 120051: iGrochalski eTail l2008t lAnnibali et all 
)08). Second, this would require that dwarf LSB galax- 
ies have an ISM that is more porous to ionizing photons 
than more massive galaxies. Dwarf galaxies with appar- 
ently porous H I distributions have been observed (e.g. HoII; 
iPuche et"aT]|1992l) . However dwarf galax ies typically contain 
a higher fraction of their mass in the ISM (Ivan Zee et all 1995b 



ISwaters et alj|2002l H06) while the lower mass densities of 
their disks suggests that that the ISM distribution should be 
"puffier". These factors should make it harder for ionizing 
photons to escap e, not easier. Indeed, Using SINGG data 
(lOey et aD F2007) argue that / esc increases with surface bright- 
ness. While we can not directly rule out a variable / esc as 
driving the fka//Fuv correlations we believe these arguments 
are sufficiently strong that a variable porosity scenario seems 
unlikely. 

5.4. Stochastic limitations 

iThilker, et aT] d2007l) and lBoissier et all d2007l) use stochas- 
tic effects to explain the often seen strong UV emis- 
sion bejwid_flie_jIatnincation radius of spiral galax- 
ies dMartin & Kennicuttl 1200 ll) . They argue that the low 
Fhq//fuv results from stars forming in events that are too low 
in mass to have been likely to form O stars. In other words, 
the IMF is normal, but the star formation is too weak to pro- 
duce the high mass stars. However, If the IMF is purely a 
statistical distribution of the masses of stars formed in a sin- 
gle event, then by averaging over many active star forming 
regions, one should recover the expected mean fka//Fuv for 
the stellar population. 

We find that the galaxies in our sample are all bright enough 
that they should have multiple O stars. We calculate this as 
follows. The minimum SFR needed to be likely to see a single 
O star is given by 



SFR m ; n — 



(to) No 



(10) 



where is the average star mass, (to) is the typical life- 

time of an O star, and N*/No is the ratio of the total number 
of stars and the stars formed above the O star mass limit for 
the chosen IMF. For a single power law IMF then 
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Where we adopt Mo = 20 .M© as the minimum mass of an O 
star. For our fiducial IMF parameters (Sec 13.21) then (M.*) = 
0.35 M Q , NJN = 1510. Table [TJ shows that for this IMF, 
and long duration star formation, the median mass contribut- 
ing to L Ha is M.+ = 54 Mq. We ado pt the main sequenc e 
lifetime of such a star as (to) = 4.0 Myr (Bressan et al.ll 19931) . 
hence the SFR must exceed SFR m ; n = 1.3 x 10" 4 A^oyear" 1 
for it to be likely to observe one O star. This ag rees very 
well with the simulations of Thil ker. et al.l d2007l) showing 
that the "stochastic limit" for observing Ha emission is at 
~ 10 -4 M-q year -1 for the same IMF parameters. 

For our fiducial stellar population, we calculate that SFR m j n 
corresponds to Mfuv = -8.29 ABmag. This is a conservative 
estimate since shorter Af would result in a fainter Mfuv. All 
the galaxies in our sample are significantly brighter than this, 
and we calculate that they should have No ^16 for continu- 
ous star formation with the fiducial IMF. However, we note 
that HIPASS J 132 1-31 has a pre-dust extinction correction 
absolute magnitude of M FUV = -9.61, the faintest in our sam- 
ple. This implies that we are only likely to see three or four 
O stars in this galaxy. This is the only source that looks like 
a complete Ha non-detection in our sample. With this few O 
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FIG. 8. — Tracks of -Fhq//fuv versus Ehc* (panel a) and (panel b) made by the burst and gasp SFHs shown in Fig.[7jusing the same convention for line 
color to distinguish event duration and line thickness to represent event amplitude. The data is the same as in Fig. [5] but displayed in gray so the model lines show 
up clearly. The model trajectories converge on the pre-event values, for which we have adopted Ehu and E« zeropoints to be consistent with typical galaxies 
having the equilibrium fHa//FUV- 



stars likely to be present at any one time, it could be that there 
is no Ha because the dust absorption correction is overesti- 
mated and we caught the galaxy at "a bad time". Otherwise, 
our target galaxies are well beyond the limit where statistical 
effects will bias /*Ha//Fuv- 

Stochastic effects may be even more important once 
plausible physical constraint s are put on the IMF. Such 
an approach is adopted b y iKroupa & Weidnerl (I2003I) and 
Weidn er & Krou pa (2005]) who consider the clustered na- 
ture of star formation, and impose the requirement that M. u 
of a given cluster is statistically allowed for the mass of 
stars formed in the cluster -M c iuster- We refer to this as 
the cluster mass constrained IMF scenario. For an underly- 
ing 7 = -2.35 (for the highest mass stars in a broken power 
law IMF following Kroupa 2001, 2002) this approach re- 
sults in a strong correlation between M. u and .M c i us ter for low 
masses (M. u oc -Mciustei- ' 67 ) that asymptotes to the true stel- 
lar M u ~ 150 Mp, for cluster s having 7W c i ustei > 6800 M Q 
(Pflamm-Altenburg et al. 2007). For a power-law cluster mass 
function A^ciustei oc Ai~ 2 the resultant IG IMF slope steepens 
to 7 ~ -2.8 for the highest ma ss stars (IKroupa & Weidnerl 
120031: IWeidner & Kroupa! 120051) . The cluster mass con- 
strained IMF scenario predicts that traditional Ln a to SFR 
conversion factors underestimate the true SFR, especially 
for SFR < 10~ 2 M.q year" 1 , up to two orders of magnitude 
higher than our estimate of SFR m ; n . The scenario can also 
account for the correlation between the luminosity of the 
bri ghtest young st a r clus ter and the SFR obse rved in galax- 



tes 



Wei dner et al.l (120041) . In a similar vein, iThilker. et al.l 



J2007F also simulate the effects of requiring a certain mass 
of stars to form in each cluster before high mass stars can 
form, and find that such a prescription could result in an or- 
der of magnitude increase in the stochastic limit SFR m ; n . This 
could provide a reasonable explanation for the low i*Ha/./Fuv 
values observed in LSB galaxies and outer disks. Indeed, 
Pflam m-Altenburg & Kroupa d2008l) show the cluster mass 
constrained IMF scenario can produce Ha edges and the non- 
linear relationship between the projected ISM density and 
Ehq: observed within galaxies. Further implications of this 
scenario for our results are discussed in Sec. 16.31 
While the cluster mass constrained IMF scenario is based 



on an IMF that has an invariant form, the resultant IGIMF is 
variable and only approaches this form for very high SFRs. 
Thus this can not be considered a constant IMF, but rather an 
explanation of what could cause the IGIMF to vary. Further- 
more, while the scenario is appealing, the physical basis for it 
is somewhat elusive. It in essence it is replacing a probabilis- 
tic interpretation of the IMF "the formation of a high mass 
stars in a low mass cluster is unlikely" with a deterministic 
interpretation "the formation of a high mass stars in a low 
mass cluster is impossible". The implication is that M. u is set 
by the available gas supply. However, the formalism is based 
on the final stellar mass of the cluster, not the gas mass of 
the progenitor. Forming or "embedded" clusters are found i n 
molecular clouds having M > 100 Mq dLada & Ladall2003h . 
sufficient to form at least one star with our adopted M. u , so 
it is not the ISM supply that is limiting M. u . In reality the 
molecular clouds fragment, and especially for low mass clus- 
ters, the star format ion efficiency is lo w and the clusters do 
not remain bound dLada & La da 2003). So in effect, nature 
may produce stars in a way that resembles the cluster mass 
constrained IMF. In Sec. 16.21 we discuss an alternative sce- 
nario for a variable IMF based on the physics of high mass 
star formation. 

5.5. Metallicity 

Metallicity affects /*h q //fuv via line blanketing, which 
affects the stellar atmosphere temperatures. We calculated 
Starburst99 stellar population models with the same param- 
eters as our fiducial population model but with metallicities 
of Z = 0.004 and 0.05 (1/5 and 2.5 times solar metallicity). 
These metallicities were used for both the isochrones and the 
spectra used to calculate the fluxes. These models result in 
Fhq//fuv = 14.0 and 8.2 respectively. Decreasing metallic- 
ity increases /*h q //fuv as the hotter stellar atmospheres re- 
sults in relatively more ionizing photons compared to the non- 
ionizing UV continuum. The models cover from SMC to su- 
per solar metallicities and result in a ~ 0.25 dex change in 
Fhq,//fuv, far smaller than the range of the /*h q //fuv varia- 
tions. Since LSB dwarf galaxies typically have lower metal- 
licities than normal galaxies, any trends with £ should act 
to reduce the observed correlations, rather than contribute to 
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FIG. 9. — Distribution of dynamical times of the sample galaxies consisting 
of single ELGs having a/b > 1.4. The dotted line shows the median of the 
distribution. 

them. Therefore, we conclude that metallicity does not con- 
tribute to the observed correlations of Fhq//fuv with surface 
brightness. 

These results hold for equilibrium models with our fidu- 
cial IMF parameters. Further exploration of models having 
different parameters show that for very steep 7 < -3.3 and 
Z = 0.004 Ffta/fpuv becomes smaller than for solar metallic- 
ity models. This is shown in Fig. [Tob discussed below. As 
metallicity decreases lower mass stars make a larger contri- 
bution to the FUV emission, and for very steep 7 values this 
becomes more important than the increased output of ionizing 
photons. 

5.6. IMF variations 

Figure [10] show results of Starburst99 models designed to 
probe the sensitivity of Fhq//fuv to the IMF parameters M. u 
and 7. The models have constant SFR over duration Af up to 
1 Gyr, and (mostly) solar metallicity. We hold other parame- 
ters fixed to the values of our fiducial stellar population model. 
In panel a. A4 U is varied and 7 = -2.35 is held fixed, while in 
panel b. 7 is varied and M. u = 100 Alois held fixed. In this 
panel we also explore cases with Z = 0.004 and Z = 0.05 and 
Af =1 Gyr in panel. Panel (c) shows the observed distribution 
of Fua/fpuY for comparison. We now consider whether the 
observed range of Fu a /fFuv values can be modeled with long 
duration star formation by varying A4 U or 7. The full sample 
(excluding upper limits) spans Fn a /fp\jv = 1.5 to 30A. Since 
the extremes may be affected by measurement errors, SFH 
variations, and stochastic effects we also consider the tenth to 
ninetieth percentile range of Fu a / fpvv from 3 to 19 A. 

For a fixed (Salpeter) 7 = -2.35 it is possible to reach the 
lowest Fhq//fuv with Ai„ ~ 20 M.q, and the tenth percentile 
with M. u = 30^0. With equilibrium star formation models 
(Af > 300 Myr) it is not possible to reach either the nineti- 
eth percentile or maximum Fhc*//fuv values. The maximum 
M u = 120 Mq allowed by the Starburst99 models results 
in an equilibrium Fn a /fyuv = 12. 9A. The curves in Fig. [TOh 
appear to be asymptotic towards high M. u so it is unlikely 
that exploring higher Ai u would help. Lowering metallic- 
ity to Z = 0.004 brings the equilibrium Fh q //fuv = 14.0 for 
M u = IOOMq (Sec. 15.51 ) helping to alleviate the problem. 



To achieve the highest observed Fh q //fuv values with solar 
metallicity and a Salpeter 7 requires Af < 10 M yr, sh orter 
than the typical dynamical times of galaxies (Sec. 15.21 ). and 
thus rarely to be found in nature. 

Varying 7 can explain a much larger range of the observed 
Fn a /fF\jv values. Equilibrium star formation with 7 = -3.5 
can account for the lowest Fhq//fuv values, while the tenth 
percentile corresponds to 7 = -3.3. The ninetieth percentile 
corresponds to 7 = -1.9, while the maximum Fhq//fuv is 
slightly larger than the equilibrium Fna/fvuy = 26 A for the 
shallowest 7 = —1.3 we considered. For Z = 0.004 the maxi- 
mum Fh q //fuv corresponds to 7 w -1 .7. 

Having ruled out other possible explanations, we conclude 
that most of the observed range of Fhq.//fuv values can most 
plausibly be explained as arising from long duration star for- 
mation with an IMF that varies at the upper end from galaxy 
to galaxy. The lowest observed values can occur by lowering 
M. u to ~ 30 M.q, steepening 7 to ~ -3.3 or an intermedi- 
ate combination of the two. It is harder to model the absolute 
highest observed Fhq//fuv values with equilibrium star for- 
mation at solar metallicity; lower metallicity can alleviate this 
problem. 

6. DISCUSSION 
6.1. Previous work 

Ours is not the first study to probe the nature of the IMF us- 
ing Ha and UV observations, nor the first to claim e vidence of 
IMF variations at the upper end. iBuat. et all (fl987h used UV 
observations from the SCAP2000 experiment to explore the 
nature of the Ha to UV flux ratio in a sample of bright spiral 
and irregular galaxies. They found a correlation in the ratio 
with morphological type, consistent to what we find: early 
type spirals have higher Fu a /fFuv compared to late type spi- 
rals and irregulars. Technical issues (multiple instruments, 
uncertain dust corrections) limited them from making strong 
claims about the IMF. They also noted a varying / esc could 
also explain their observed trend. 

Recently, IHoversten & Glazebrook (2008) deduced varia- 
tions of the upper end of the IMF using a very large (~ 140K 
galaxy) sample from the Sloan Digital Sky Survey (SDSS) 
and employing an Ha equivalent width ver sus optical color di - 
agram, following the method pioneered by Kenn icuttl d!983l) . 
They find that low luminosity galaxies have much lower 
EW(Ha) than expected for their colors and argue that this is 
due to the upper end of the IMF varying systematically with 
galaxy mass. The variations are not monotonic with lumi- 
nosity, but nevertheless similar to our results: low luminosity 
galaxies have less massive stars than higher luminosity galax- 
ies. Their result is solid, but since optical bands are involved, 
much lower mass stars and longer timescales were required in 
their modeling to rule out other explanations. 

Our contribution t o the field is to return to the basic method 
of IBuat. et al.l ([1987) with a larger sample and consistent dust 
absorption corrections. By using Ha and FUV measurements 
we isolate the upper end of the IMF, and hence our constraints 
are more direct than Ha - optical color comparisons. A new 
insight from our work is that the surface brightness, or more 
likely surface mass density, seems to be more important for 
driving the IMF variations than luminosi ty or mass as sug- 
gested by IHoversten & Glazebrookl (l2008h . Our use of use of 
H I selection and measurement of Ha and UV images results 
in an extended range of optical surface brightness compared to 
SDSS (M06), allowing this finding, while our inclusive sam- 
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FIG. 10. — Starburst 99 model predictions of Fho/Zfuv for stellar populations forming at a constant star formation rate are shown in panels (a) and (b). Panel 
(c) shows the distribution of observed FH a //FUV values for comparison. In panel (a) the IMF slope 7 is held constant at the Salpeter (1955) value while M. u is 
varied. In panel (b) the upper mass limit is held constant at A4 U = 100 Mq while 7 is varied. In these panels the dots indicate the _M„ and 7 values used in the 
solar-metallicity (Z = 0.02) models, while solid lines connect the models for star formation durations of 0, 10 Myr, 100 Myr, and 1 Gyr (from top to bottom). The 
dashed (purple) and dot-dashed (green) lines in panel (b) indicate tracks with a metallicity Z = 0.004 and 0.05 respectively and a duration of 1 Gyr. The broken 
horizontal line in all panels shows the expected Fhc<//fuv for our fiducial stellar population model (eq.f5J. 



pie selection eliminates the possibility that our results stem 
from a non-representative sample. 

6.2. A scenario for the observed correlations 

The belief that the form of the IM F is universal i s largely 
based on observations of star clusters. iRroupal d2001l) showed 
that apparent variations in 7 observed in young star clusters 
could result purely from stochastic effects due to the limited 
number of massive stars. A universal IMF is then a logi- 
cal conclusion if all stars form in star clusters as argued by 
lLada&Ladald2003l) . 

Recent work suggests that this paradigm should be re- 
examined. Theory and simulations indicate that the highest 
mass stars (> 10 A4q) require a dense cluster environment to 
form efficiently. It is hard to form these stars by simple col- 
lapse and fragmentation - the stars "turn on" and stop form- 
ing before high masses can be assembled. One way around 
this is by the process of competitive accretion whereby young 
proto-stars "steal" material from each others envelopes as well 
as accreting it from the dense ISM at the bottom of a cluster's 
potential well (iBonnell. et al.l2 004). This allows rapid assem- 
bly of massive stars but requires very dense env ironments, like 
star c lusters, and especially their dense centers (IBonnell. et al.l 
2004). While all stars may form in clusters, not all clusters 
are alike. Clusters with a wid e range of densitie s are known 
many of which are unbound dLada & La da 2003). Ephemeral 
unbound clusters are not likely to maintain a high enough den- 
sity for long enough to have the collisions necessary to form 
the most massive stars. 

Pressure p lays a key role in re g ulating the phases 
of the ISM dMcKee & Ostrikerl [l977t IWolfire et al.l 120031; 
iBlitz & Ro solowskv 200g), and is the most likely phys- 
ical property for driving the observed Fhq//fuv correla- 
tions. While the neutral phases (cold and warm) are in pres- 
sure equilibrium with each other, the molecular ISM is self 
gravitating and supported by turbulen ce from star formation 
feedba ck dDopita & Sutherland! 120031) . IBlitz & Rosolowskv! 
(2006) show that the ratio of molecular to neutral ISM in 
galaxies has a nearly linear relationship with the mid-plane 
pressure, f m id, expected from hydrostatic equilibrium of a thin 
gas layer embedded in stellar disk with a significantly larger 
scale height. This leads to the expectation that 



(13) 

T, g is the ISM 

mass density, a g is the ISM velocity dispersion, and h+ is the 



where is the stellar mass density in the disk, ^ g 



stellar scale height. The terms in parenthesis are not expected 
to vary much within galaxies nor from galaxy to galaxy. 
varies strongly between galaxies and to first order is oc 'Er. 
The molecular to neutral mass ratio in the ISM, and therefore 
the fraction of the ISM available for star formation, should 
be highly dependent on £«, because it effectively is tracing 

p ressure. 

lElmegreen & Efremovl d 19971) and lElmegreenl d2008bl) 
show that tight bound clusters, and hence O stars, should 
preferentially form in a high pressure environment while at 
lower pressures unbound and loose clusters form (with less 
O stars). While the pressure they are referring to is that in- 
ternal to the molecular clouds, P m id provides a floor to this 
pressure and t hus should have some beari ng on the final inter- 
nal pressures (Dopita & Sutherland 2003). By this reasoning, 
£r traces P m ;d which determines the likelihood that clusters 
are formed bound, thus regulating the O/B ratio and therefore 

•Fhq//fUV- 

The link between compact cluster formation and 
surface brightnes s is well established observationally. 
lLarsen & Richflerl (|2000) and iLarsenl d2004l) used ground 
based and HST images of spiral galaxies to show that the 
fraction of U band light in star clusters corr elates with 
the su rface brightness of the host. Similarly, Bill et et aT] 
(2002) show that the fraction of U band light in star cl usters 
scales almost linearly with Ssfr- Me urer et alj d 1995b used 
Hubble Space Telescope NUV images to show that there is 
a correlation between fraction of NUV light in the form of 
star clusters and underlying surface brightness in starburst 
galaxie s . Whi le the correlation is weaker than that shown by 
ILarsenl d2004l) . the surface brightness range is smaller, with 
the starbursts corresponding to the high intensity, high cluster 
fraction extension of the sequence of Larsen. 

Internal variations in fka//Fuv within a galaxy may also 
be related to pressure variations. For e xample, M83 shows 
a sharp decline in its i*k Q //Fuv ratio dThilker et alJ 12005b 
corresponding to the Ha truncation radius identified by 
iMartin & Kennicuttl d2001l) . The FUV surface brightness pro- 
file continues well beyond this radius, with no apparent sign 
of a truncation . Simi lar results are found in a larger sample by 
iBoissier et alj d2007[) showing that H II edges, generally are 
not seen in the FUV. Martin & Kennicuttl (l2001h show that the 
H II edges corresponds well to where S s drops below the crit- 
ical d ensity needed for disk self-gravity (M artin & Kennicuttl 
1200 lb . When a disk becomes gravitationally unstable it first 
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excites low order modes: bars and spiral arms. These will in- 
duce spiral shocks which will increase the ISM pressure and 
thus facilitate the formation of star clusters and the efficient 
formation of O stars. Thus star-formation edges may be a 
byproduct of the sharp-pressure increase associated with spi- 
ral arms in massive unstable disks. 

Our scenario is somewhat complementary to the more 
mathematically developed cluster mass constrained IMF sce- 
nario of Kroupa and Weidner discussed in Sec. 15.41 Indeed, 
they could be compatible to the extent that our scenario offers 
a physical basis for small unbound clusters having a truncated 
IMF. 

6.3. Implications of a variable top end of the IMF 

If the upper end of the IMF is not uniform then there are 
several important astrophysical implications. First, the star 
formation rate (SFR) estimated for a galaxy depends on the 
tracer used. While Ho and UV based SFR may agree for the 
bright galaxies, Ha based SFRs will underestimate the contri- 
bution of LSB galaxies to the cosmic star formation density of 
the universe. Even if an alternate explanation of the Fh q //fuv 
variations were found (e.g. / esc variations), SFR estimates will 
still be under-estimated from Ha fluxes in LSB galaxies. Op- 
tically selected surveys of star formation will have results that 
depend on their surface brightness selection limits. Perhaps 
little notice has been paid to the problem of low Fn a /fpuy 
values because most star formation surveys have concentrated 
on normal and HSB galaxies where the ratio is closer to nor- 
mal. UV based SFR estimates should be more secure than Ha 
estimates. However, there is no guarantee that IMF variations 
stop at the O - B star divide. For the cluster m ass constrained 
IMF scenario Pflamm-Altenburg et al. (2007) show that when 
the SFR is very low (SFR < 10" 3 M yr _1 ) that the IGIMF de- 
viates from the canonical IMF down to M.+ of just a few A4q 
and that standard conversion factors between L# Q and SFR 
could underestimate the SFR by orders of magnitude in such 
cases. However, at the low mass end, below ~ 0.8 the 
IMF appears to be constant (Feltz ing et al.l ll999: Wys e et all 
2002). More work needs to be done to determine where be- 
tween the masses of O stars and the sub-solar range that the 
IMF switches from being variable to constant. 

Second, a varying IMF also implies that the feedback 
of the young stellar populations onto the ISM in terms of 
energy output and the return of chemically enriched ma- 
terial will also vary, with LSB galaxies having less feed- 
back of energy and mass per mass in stars formed. The 
metallicity of the returned ISM will also be affected. LSB 
galaxies will have a lower [O/Fe] ratio and higher [N/O] 
ratio than normal or HSB galaxies. While such abun- 
dance anomalies are ob served in LSB dwarf galaxies (e.g. 
iTautvaisiene et alJ l2007h the abundance patterns are usu- 
ally explained by galactic winds that preferentially expel 
oxygen rich ejecta in galaxie s with weak gravitational po- 
tentials dPilvuginl 1 1 992l 1 1 993t iDe Young & Heckmanl ll994t 
iMac Low & Ferraralll999l: lD , Ercole & Brighentilll999l) . Un- 
doubtedly such winds occur (e.g. NGC 1569: Waller 1991; 
NGC 1705: Meurer et al. 1992), but they would be harder 
to produce in a LSB dwarf compared to a HSB blue com- 
pact dwarf. Instead of removing excess metal content in a 
galactic wind, LSB galaxies m ay just not produce the metals. 
Recently, iKoppen et alJ (|2007) were able to model the func- 
tional form of the mas s-metallicity and mass yield relations of 
dTremonti et alj[2004b using a cluster mass constrained IMF 
model described above (Sec 15.41 ) effectively demonstrating 



this scenario is plausible. 

Measurements of the SFH of nearby galaxies from color 
magnitude diagrams (CMDs) can also be misinterpreted. In 
such an analysis, one typically assumes that stars form with 
a Salpeter IMF fully populated up to ~ 100M S . If M u 
is lower, or 7 steeper than assumed, then one could spu- 
riously infer recently truncated or declining SFHs. Indeed 
many galaxies in the nearby universe show very few field 
high mass stars . These include LSB dwarf irregu lars (e.g. 
iTosi et al.l fl99lt iGreggio et alj|1993t iTolstovll 19961). but also 
some blue compact dwarfs (IGreggio et al.ll998HCannon et alJ 
120031) . The CMP results can then be interpreted as a SFH 
that is "gasping" ( Marconi et al. 1 19951) or being in a "post- 
burst" phase dCannon et al.ll2003 ). Gasping could also result 
from a bias inherent in the star formation scenario: stars are 
harder to identify in crowded clusters, hence there is a bias 
against finding stars with life times shor ter than the dissolu- 
tion time of clusters (Marconi et al. 1995). However, there are 
also cases where the gasping deduc ed from CMDs is consis- 
tent w ith weak Ha measurements (lTolstoyi ri996: Co le et all 
120071) . Claims based on CMD analysis that star formation is 
"contracting" or pro gressively halting fro m large radius to an 
inner starburst (e.g. Skill man et alJ l2003) may also be cast in 
doubt if there is a radial variation of the IMF. 

7. SUMMARY 

We have looked at correlations of the integrated Ha to 
far UV flux ratio fka//Fuv with other global parameters in 
a sample of galaxies that well represents the full range of 
star formation properties of galaxies in the local universe. 
There are strong correlations of Fh q //fuv with optical surface 
brightness in both Ha and the R band. Weaker but significant 
correlations are found with R band luminosity, rotational am- 
plitude as well as morphology. Thus the systematic variations 
of Fua/fmv are part of the family of observational galaxy 
scaling relations, such as the Hubble sequence, the Tully- 
Fisher relation, the mass-metallicity relation, the luminosity- 
surface brightness relation and the star formation law. How- 
ever, unlike other scaling relations, the ^ka/ZFUv correlations 
appear to be more closely related to high mass star formation 
intensity or stellar mass density than halo mass. 

We examined a variety of plausible explanations for the root 
causes of these correlations and ruled most of them out as fol- 
lows. 1. While dust correction decreases Fhq/Zfuv, the effect 
is not large enough to have spuriously created the correlations, 
and can not fix the problem of low fk a //Fuv values which 
exists even before dust correction. 2. Recent variations in the 
SFR would need to have had large amplitudes (factor > 10) 
and short duration (< 100 Myr) to account for the range of 
Fhq//fuv- The normal undisturbed morphology of the major- 
ity of our sample belie such extreme events in their recent his- 
tories. Furthermore, variations in the SFR can not account for 
the correlations of Fhq//fuv with both Sh q and Invok- 
ing SFH to explain the low Fu a /fFuv in LSB galaxies would 
require all LSB galaxies to be experiencing coordinated gasps 
in their SFR which violates the Copernican principle. 3. Low 
Fn a /fpuv values can occur if ionizing photons escape. How- 
ever, one would then expect to see naked O stars in Hubble 
Space Telescope images of the nearest LSB galaxies, and typ- 
ically they are not found. Furthermore, LSB galaxies tend to 
have a higher gas fraction, and should have thicker disks than 
normal galaxies which argues against them having a higher 
/esc- 4. Stochastic effects can not cause the observed trends 
since we are measuring over entire galaxies, and they all are 
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sufficiently luminous in the UV, even without dust correction, 
that they should contain multiple O stars. The related cluster 
mass constrained IMF scenario of iKroupa & Weidnerl (2003) 
and Weidn e~& Kroupal (|2005) may explain the results, but in 
effect it does not describe a constant IMF but a variable one. 
5. Metallicity has only a fairly small effect on Fn a / /ftjv with 
lower metallicity acting to increase Fn a /fpu\ for a normal 
IMF. This is also in the wrong sense since LSB galaxies tend 
to have low metallicity and low ^hq/Zfuv- 

This leaves the most plausible explanation of the results 
to be a varying upper end of the IMF. The observed range 
Fhq//fuv can be explained if the upper mass limit M. u varies 
between 30 and 120 M.q, or the IMF slope 7 varies between 
-1.3 and -3.3. There is not enough information to decide be- 
tween these options or rule out variations in both parameters. 

We sketch a scenario that can explain the observed corre- 
lations as arising from the pressure regulated formation of 
tightly bound star clusters. These are the sites of the formation 
of the highest mass stars. A steeper or truncated IMF occurs 
when stars form in loose or unbound clusters. Regions of high 
pressure favor the formation of bound clusters resulting in an 
increased O/B star ratio and thus higher Fhq//fuv- R band 
surface brightness traces pressures since it is the established 
stellar populations that dominate the disk potential and hence 
determines the mid-plane pressure at hydrostatic equilibrium. 

There are several major implications of these results. The 
SFR estimated from Ha and FUV measurements are not nec- 
essarily the same. This holds even if the correlations result 
from a variable / esc instead of IMF variations. If the IMF is 
to blame, then the return of energy, momentum, and metals to 
the ISM from high mass stars is decreased in LSB galaxies. 
This provides an alternative driver for the mass metallicity 
relation other than galactic winds. Allowing IMF variations 
makes the interpretations of Color-Magnitude Diagrams dif- 
ficult: an observed lack of high mass stars may be due to a 
truncated IMF rather than being due to a recent decrease in 



the star formation rate. Hence interpretation of CMD results 
showing gasping/post-burst SFHs or radial truncations of star 
formation are cast in to doubt. 

Further work is required to pin-down the nature of the 
Fhq//fuv variations. We are currently examining internal 
Fn a /fF\jv variations within our sample galaxies in order to 
determine the extent to which the IMF variations are local. 
More observations of LSB galaxies are needed to firmly rule 
out a variable / esc as the cause of the variations. If LSB galax- 
ies do have a high / esc this would imply that either their H II 
regions no longer follow case B ionization, or that they have 
many naked O stars. Spectroscopic tests for either scenario 
should be possible. Careful CMD analyses, including UV 
measurements, of HSB and LSB galaxies can be used to deter- 
mine how far down the IMF the variations may occur. Com- 
parison of abundance ratios with surface brightness may pro- 
vide constraints on the relative importance of a truncated IMF 
and starburst winds for determining the abundance patterns in 
dwarf galaxies. 
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